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Abstract. This study starts with balances deduced by Baumert and Peters 
(2004, 2005) from results of stratified-shear experiments made in channels 
and wind tunnels by Itsweire (1984) and Rohr and Van Atta (1987), and 
of free-decay experiments in a resting stratified tank by Dickey and Mellor 
(1980). Using a modification of Canuto's (2002) ideas on turbulence and 
waves, these balances are merged with an (internal) gravity-wave energy 
balance presented for the open ocean by Gregg (1989), without mean-flow 
shear. The latter was augmented by a linear (viscous) friction term. Gregg's 
wave-energy source is interpreted on its long-wave spectral end as internal 
tides, topography, large-scale wind, and atmospheric low-pressure actions. In 
addition, internal eigen waves, generated by mean-flow shear, and the aging 
of the wave field from a virginal (linear) into a saturated state are taken 
into account. Wave packets and turbulence are treated as particles (vortices, 
packets) by ensemble kinetics so that the loss terms in all three balances have 
quadratic form. Following a proposal by Peters (2008), the mixing efficiency 
of purely wave-generated turbulence is treated as a universal constant, as 
well as the turbulent Prandtl number under neutral conditions. It is shown 
that: (i) in the wind tunnel, eigen waves are switched off, (ii) due to remotely 
generated long waves or other non-local energy sources, coexistence equilibria 
of turbulence and waves are stable even at Richardson numbers as high as 
10 3 ; (iii) the three-equation system is compatible with geophysically shielded 
settings like certain stratified laboratory flows. The agreement with a huge 
body of observations surprises. Gregg's (1989) wave-model component and 
the a.m. universal constants taken apart, the equations contain only one 
additional dimensionless parameter for the eigen-wave closure, estimated as 
Y « 1.35. 



1. Introduction 
1.1. General 

In geophysical flows, turbulence is ubiquitous. Today 
turbulent engineering flows may already be simulated 
using DNS, i.e. fully resolving scales down to the small- 
est ones, provided they are large compared with molec- 
ular cluster scales. However, for ocean and atmospheric 
flows this will remain impossible even in the foreseeable 
future. Here particularly the dominating stably strati- 



fied flows, i.e. all forms of coexistence of turbulence and 
internal waves, represent a specific observational and 
theoretical challenge. Without a better understanding 
of the fundamental physics of these processes all our re- 
gional to global models of weather and climate remain 
incomplete. 

Oceanic, atmospheric and stellar turbulence has 
challenged a number of prominent scientists and re- 
search groups since long, in the recent decade namely 
Woods (2002), Galperin et al. (2007), Canuto et al. 
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(2008), and Zilitinkevich et al. (2008), also Kantha and 
Carniel (2009). They all emphasized a major contra- 
diction between observational experience and existing 
theories: 

(i) For controlled stratified shear flows in the labora- 
tory exists undoubtly a critical gradient Richard- 
son number of R b g < 0{^/±) above which tur- 
bulence dies out (Itsweire, 1984; Itsweire et al., 
1986; Rohr and Van Atta, 1987; Rohr et al., 1987; 
Rohr et al, 1988a, b; Van Atta, 1999). Also in 
the field the qualitative and strongly (but not to- 
tally) limiting role of R g is without doubt (Pe- 
ters et al., 1988) and became specifically visible 
with the advent of Lagrangian floaters in the La- 
grangian time spectra of turbulence and wave-like 
fluctuations where Q — N marks a sharp divide: 
to the left a flat wave spectrum, to the right a 
Kolmogorov time spectrum (D'Asaro and Lien, 
2000a, b) . Also the existing theories (Richardson, 
1923; Miles, 1961; Howard, 1961; Hazel, 1972; 
Thorpe, 1973, Abarbanel et al., 1984; Baumert 
and Peters, 2004) point all into the same direc- 
tion. 

(ii) Geophysical flows are reported to exhibit not al- 
ways, but more than often significant stable tur- 
bulence levels and mixing capabilities at R g 3> 
y 4 (Peters et al., 1988; Canuto, 2002; Poulos et 
al., 2002; Nakamura and Mahrt, 2005; Grachev 
et al., 2005, 2006). Also Peters and Baumert 
(2007) report problems in validating a K-e tur- 
bulence closure against comprehensive cstuarine 
microstructure measurements in the Hudson river. 
For tidal phases with weak shear the observed tur- 
bulence levels exceeded the model values by 2 to 
3 orders of magnitude. The dynamic time lag of 
the turbulent state variables in the M 2 tide (about 
one hour, see Baumert and Radach, 1992) can be 
excluded as a source of the deviations because the 
model naturally contains this effect. 

The text below makes an attempt to resolve the 
sketched contradiction using ideas of Woods (2002 and 
literature cited therein) and Canuto (2002) regarding 
the role of waves and billow turbulence. We identify 
internal gravity waves (so far neglected in most other 
studies) as the potentially responsible phenomenon. 

For a better understanding of the language used in 
this text we first introduce notational conventions and 
then continue with an overview of major physical pro- 
cesses within the world of stratified shear turbulence 
and internal gravity waves. 



Below, we always means the author and the dear 
reader in a dialogue. 

1.2. Setup and notation 

For simplicity we focus on a simple non-trivial situation, 
a spatially one-dimensional "channel" flow with velocity 
component U in horizontal (x) direction, with variation 
of U along the vertical, z (pointing upwards), and with 
Eulerian density and icopycnal coordinate fluctuations. 

The decomposition of our flow field into mean and 
fluctuations reads as follows: 

U{z,t) = (U)+u(z,t)+u'(z,t), (1) 
W(z,t) = (W) + w(z,t) + w'(z,t). (2) 

Variables with tilde are small-scale short-wave compo- 
nent^ present only under stratified conditions. Primed 
variables denote turbulent fluctuations. Both fluctuat- 
ing components vanish in the mean. 

The turbulent kinetic energy, TKE or JC, is defined 
as follows, 

K = \(u' 2 +v' 2 + w' 2 ). (3) 



A wave's total energy, £, is the sum of potential, 
(N 2 ( 2 )/2, and kinetic energy, (u 2 +w 2 )/2 (Gill, 1982): 



£ = \{u 2 +w 2 + N 2 ~C 2 ). 



(4) 



Here u points into the direction of wave propagation so 
that in a plane wave v = 0. Average (.) is taken at least 
over one wave period which "by definition" is longer 
than the characteristic turbulent time scale. In the fol- 
lowing the dimensionless gradient Richardson number, 
R g , plays a central role. In a stratified shear flow as 
above, shear is given by 



(S) = 



dz 



(5) 



while stratification is characterized by the Brunt- Vai- 
sala, frequency squared, 



(N 2 



9_ d(p) 
(p) dz 



(6) 



(p) is the background density field. 

The gradient Richardson number, R g , characterizes 
the dimensionless ratio of the two aspects, shear and 
stratification: 



R<> 



(N 2 )/(sy< 



(7) 



1 actually wave packets 



Baumert: Turbulence by shear & waves 



3 



Below the averaging operators are mostly omitted for 
brevity of notation. But fluctuations are consequently 
labeled either by tilde (wave-like) or prime (turbulent). 

Based on R g , under controlled laboratory conditions 
where the linear eigen waves leave the experimental 
site before quadratic saturation and a feed back into 
the TKE pool can happen (wind tunnel of Van Atta) , 
the following hydrodynamic regimes for horizontally ho- 
mogeneous flows are found (Baumert and Peters, 2004, 
2005): 

(a) R g < Rg = 0: unstable and neutral stratification, 
convective turbulence, no internal waves at all. 

(b) = Rg < R g < R b g = V 4 : stable stratification, 
shear-dominated growing turbulence, coexistence 
of turbulence and internal waves. 

(c) 7 4 = R g < R g < R g = Y2: stable stratification, 
wave-dominated decaying turbulence, coexistence 
of turbulence and internal waves. 

(d) 72 = Rg < R g : stable stratification, no coexis- 
tence of turbulence and waves, waves-only regime. 

Under those conditions the turbulent Prandtl number 
cr is (Baumert and Peters, 2004, 2005) 

a -- v - 1 - (iv/op (8) 

The above values for the critical numbers R a gl R h g and 
Rg hold for the asymptotic case Re — > 00. 

The situation is different if we leave the laboratory 
wind tunnel and consider the open ocean, stratified 
rivers or the stably stratified atmosphere where turbu- 
lence is not only locally generated through local mean- 
flow shear but also through the action of space-filling 
(non-local) spectra of internal gravity waves (IGWs). 
These are generated e.g. by tidal forces, possibly at re- 
mote places, arriving at our point of interest along var- 
ious pathways. 

2. Major physical interactions in stably 
stratified oceans and atmospheres 

Fig. [^schematically presents the major interactions be- 
tween the energies of mearj^] and fluctuating motion 
component^] Fig. [2] does the same for r.m.s. vorticity 
and turbulent viscosity, v. The latter connects K, and 17 



2 Mean-flow kinetic energy, Ya (U) 2 , MKE; Internal tides 
3 TKE, K; r.m.s. vorticity, Q; WKE. 



with the mean flow through the so-called Kolmogorov- 
Prandtl relation in the following form (Baumert and 
Peters, 2004, Baumert, 2012): 

v = K/-k$1. (9) 

The TKE dissipation rate, e, transforming TKE into 
heat, is defined as (Baumert and Peters, 2004, Baumert, 
2012) 

e = JCtt/ir = vtt 2 . (10) 

The buoyancy flux, B, transforming TKE into back- 
ground potential energy, PEt,, may be expressed using 
the eddy diffusivity, /i, as follow^] 

B = -*<«/ p')=fiN 2 , (11) 

where u, is related with the eddy viscosity v through er, 
the turbulent Prandtl number function: 

fi = vja or a = v I [i . (12) 

Finally, the molecular heat flux, transforms internal 
energy into background potential energy, PEt,. 

There is further the rate II, which transforms the 
energy of internal tides and wind-generated large-scale 
internal motions over a long chain of various friction- 
poor wave-wave interactions into short-wave internal- 
wave energy, and there is the rate P, which transforms 
the energy of internal-waves spectra into TKE, mainly 
by breaking (a shot noise process), and to a low degree 
by wave shear. 

The mean-flow shear, S, controls TKE production 
through P — and internal-wave generation through 
ty. In the (conventional) neutrally stratified case, iV 2 = 
0, the total loss term in the mean-flow kinetic-energy 
balance (MKE) is — P, which is at the same time the 
only source of small-scale energy: 

P = vS 2 . (13) 

But for N 2 > small-scale energy means the sum of 
turbulence and waves, P — (P — \P) + , 5. 

The shear governs not only TKE and IGW produc- 
tion but also the generation of r.m.s. vorticity. But vor- 
ticity is also influenced by short internal-gravity waves 
through the term S: in the shearless IGW-dissipation 
case (e.g. Gregg, 1989) we have also vorticity genera- 
tion. 



4 The buoyancy flux refers to purely shear-generated tur- 
bulence (for details see Subsection |5.1.2| and Baumert and Peters, 
2004, 2005). 
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Figure 1. Interaction of major forms of fluid- 
mechanical energy in stratified oceanic, atmospheric or 
stellar shear flows. 'J' is the flux of eigen-wave energy 
mentioned in the text. II is the energy flow which corre- 
sponds to remotely generated waves. The wind tunnel 
allows to cut off the flux P from the the wave-energy to 
the TKE pool because the wave spectrum cannot reach 
a saturated state. In the long-term equilibrium holds 
P = (1 - + II) where / is the fraction of direct 
wave dissipation and tends to zero if the wave spectrum 
approaches saturation. / may become relevant for short 
wave ages. 
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Figure 2. Interaction of major vorticity-controlling 
motion components in stratified oceanic and atmo- 
spheric shear flows. The r.m.s. vorticity fi is governed 
by the mean-flow shear, S, and a wave-induced pseudo 
shear, 5. The latter is controlled by eigen waves (Y) 
and remotely generated long waves (X). For details see 
text. 



we mostly mean for brevity TKE (JC), i.e. turbulent ki- 
netic energy. If we talk here about vorticity we similarly 
mean for brevity the r.m.s. turbulent vorticity. For its 
detailed mechanical interpretation we refer to Baumert 
(2012). 

The precise meaning of £1 is actually the module of 
the vorticity frequency. The module of the vorticity 
itself, u>, is to = 2 it il. The enstrophy is w 2 /2. 

3. Balances for small-scale motions 

Before we start to write down transport equations for 
primary variables we first discuss their nature and their 
balances. For this aim we look again at Fig. [T] and 
there namely on the two boxes in the middle row with 
the names TKE and Internal waves. 

The left box, TKE, is fed by two components: by 
mean-flow shear in the form of P — ty, and by internal- 



wave breaking, P, see ( 15 ) below. But it generates heat 
by e via ( 10 ), exports buoyancy B via ( 11 and gener- 
ates by ^ in the course of aging and subsequent collapse 
short internal waves, when its time scale approaches the 
internal- wave period (Baumert and Peters, 2005). 

Turbulence collapsing into waves is an event mostly 
bound to certain special conditions (e.g. Dickey and 
Mellor, 1980; D'Asaro and Lien, 2000a, b). Under 
smooth and almost-equilibrium conditions the TKE 
may collapse and generate waves which then after aging 
saturate and feed their energy back into the TKE pool, 
as part of P. In the following ^ is therefore mostly 
neglected. 

The right box, Internal waves, is fed by three com- 
ponents: by large-scale wave sources, II, by the rate of 
eigen-wave generation, 4", and by the collapse rate dis- 
cussed already above and neglected further below. This 
box further exhibits two relevant losses: the linear wave 
friction, L, 

l = c 1 £, (14) 



and the quadratic wave friction, P, 
P = c 2 £ 2 . 



(15) 



The latter is mainly caused by wave breaking and domi- 
nates the wave-energy balance, but not completely. The 
loss ternj^ ( 14 ) and ( 15 ) follows arguments developed by 
Gregg (1989) which we augment as follows. 

From a purely theoretical point of view the phe- 
nomenon of wave breaking tells us nothing about its 



When we use the word energy in the present context, 



5 generating thus background pote ntial energy, PE^,. 
6 see also the wave-energy balance disk below 
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kinetics. It is the processes before breaking occurs 
which form the bottleneck. If we accept the idea that 
not only vortices (actually: vortex-dipole filaments, see 
Baumert, 2012) are particles which move in space un- 
til collision, then the interpretation of wave-energy dy- 
namics in terms of particle dynamics does not come as 
a surprise, in particular in view of the billow turbulence 
discussed by Woods (2002). 

Due to the general particle- wave dualism of field the- 
ories, which is a well established concept in classical 
continuum mechanic^ also wave packets can be treated 
as particles moving with their group velocity until colli- 
sion. After collision they either move ahead, or change 
their paths, or, with low probability, their energy is dis- 
sipated by breaking in dissipative patches (billows). In 
contrast to vortex kinetics, here the probabilities are 
not symmetric. In both cases (vortex dipoles and wave 
packets) the collision events are highly intermittent. 

Besides its phenomenological basis, the quadratic 



term in (15) might therefore have deeper roots or at 



least analogies in the kinetics of reactive particle "gases" 
or molecular reactions in fluids where particle-collision 
probabilities follow the product of their spatial densi- 
ties. For collisions between particles of same kind, the 
quadratic collision term is thus a logical consequence. 

Fig. [2] shows the major feedback loop between the 
small-scale motions and the mean flow which eventu- 
ally smoothes the flow through the turbulent viscosity. 
The central left box, r.m.s. vorticity, is fed by two 
components: by the mean-flow shear, S, and by the 
internal-wave field through the pseudo shear, S. The 
vorticity f2 itself controls together with the TKE JC the 
turbulent viscosity via ([£]) which eventually smoothes 
the flow. 

The spectral signatures of shear-generated fluctua- 
tions and wave-wave interactions differ qualitatively. 
While shear influences vorticity directly by prescrib- 
ing a time scale cx S^ 1 , the long- wave sources (IT) of 
IGW spectra do it more indirectly via a longer chain of 
wave-wave interactions cascading down to critical fre- 
quencies around N. This implies that each of the two 
mechanisms needs "his" closure. 

For conditions of homogeneou^] shear, stratification 
and wave fields the above three balances can be so far 
formulated as follows (compare with Figs. [I] and [2]): 

dn i fs 2 ~ 9 \ .... 

(16) 



dt 



7 It became most famous in quantum mechanics and has been 
somewhat monopolized there. 

8 in horizontal and vertical direction 



dJC 

~dt 
d£ 

~dt 



{P-ty-B-s + P, 

(n + *) -L-P. 



(17) 
(18) 



With the exception of II, all variables in (16 17 18) 
have purely local character. 



In (17) we used the relation 

B + * = vN 2 /(T Q 



(19) 



derived and discussed in greater detail by Baumert and 
Peters (2004, 2005). Here er = 1 / 2 is the value of the 
turbulent Prandtl number function a for the case of 
neutral stratification (R g = or fi — > oo). With |8]) we 
have 



00 



ft J 



(20) 



Note that the system ( 16 - 18 1 is not one of the com- 



mon three-equation models used in traditional turbu- 
lence modeling. The focus of this modeling branch is 
directed on higher and higher orders of the closure equa- 
tions, e.g. equations for second and third moments and 
so forth, all derived from the Navier-Stokes equation. 



In later sections we will see that the system ( 16 



18) is dynamically stiff. Under spatially homogeneous 
conditions well-defined steady-state solutions exist, but 
perturbations of this state are connected with very dif- 
ferent characteristic relaxation times of the systems 
components. In particular, f2 is the component relaxing 
fastest into what we call structural equilibrium. The 
TKE (JC) relaxes significantly slower into a new state 
while the wave energy, £ , relaxes extremely slowly. The 
wave-energy pool needs a longer spin-up time because 
the flux of shear-generated wave energy, 'J, consists of 
random linear wave packets which simply need time "to 
meet and break" , of course by chance. 

For later use we introduce here the "viscous fraction" 
/ of the total energy loss of the wave pool towards the 
heat pool, 

f=^ (21) 



L + P 



with the obvious property 



(22) 



The fraction 1 — / describes a "wave age" , i.e. the rel- 
ative wave-energy loss into the TKE pool. Clearly, 



0</<l. 



(23) 
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Figure 3. Linear wave packet, hypothetically gener- 
ated by shear. Breaking occurs by chance superposi- 
tions and results in dissipative patches or billows. 




Figure 4. Solid: laminar (linear) loss fraction of wave 
energy in the course of saturation, f(t) — L(t)/(L(t) + 
P(t)); dashed: "wave age", 1 — f(t). 



We note in passing that the steady state is not the 
only dynamically invariant state of TKE and waves. 
Also the state of exponential evolution (Van Atta, 1999) 
in wind tunnels belongs to this class. This state, taken 
as a reference, has also the property that perturbations 
relax back towards reference. 

In our present situation where we deal with relax- 
ation times orders of magnitudes apart, the use of the 
so-called Tikhonov principle seems to be helpful. It 
means to concentrate on processes with moderate relax- 



ation times. ( 16 ) is so fast that it can be taken as being 



always in structural equilibrium. ( 18 1 is so slow that 



its time derivative is small compared with the source 
and sink terms at the right-hand side and can be ne- 
glected. We are thus left with only one differential and 
two algebraic equations. However, in the case of very 
stiff algebraic equations it is sometimes more useful to 
apply the method of non-stationary embedding. Here 



it would mean to re-establish the character of ( 16 1 and 



(18) as differential equations and to seek the stationary 



4. Special cases 

Below we discuss some special cases of our general sys- 
tem) 16 -18): the neutrally stratified case (N = 0), the 
stratified (N 2 > 0) but geophysically shielded (n = 0) 
case, and the stably stratified wind tunnel. 

4.1. Neutral stratification, N = 

Homogeneous shear means constant shear along the 
horizontal and vertical axes. The assumption N = 
means neutral stratification. Internal gravity waves of 
any kind are not supported by the fluid so that all terms 



in ( 18 1 vanish, together with this equation. Further, the 
terms ^, U, P, S are zero such that eventually (16 
171 look as follows: 



dn 

~dt 
dJC 
~dt 



1 

7T 
P 



s 2 



- n- 



(24) 
(25) 



These equations correspond to a clear mechanistic in- 
terpretation of turbulence as dipole chaos in the sense 
of a two-fluid approach (excitons in form of quasi-rigid 
vortex tubes made of inviscid fluid, and the materi- 
ally identical inviscid but not excited fluid between 
the tubes) derived by Baumert (2005 - 2012). As a 
byproduct, this theory gives von Karman's number as 
k= 1/V2tt = 0.399. 

4.2. Stable stratification, N 2 > 0, n = 

Here we mean stratified shear flows under idealized lab- 
oratory conditions where the role of tides and external 
geophysical influences are excluded, i.e. n = 0, which 
implies also S = 0, to be discussed later. These condi- 
tions are called below "geophysically shielded" . 

Homogeneity means here constant N 2 and S 2 in the 
horizontal plane and on the vertical axis, which is not 
easy to realize in a laboratory. But on a simplistic the- 
oretical level we can get some insight when we consider 
only the stationary system (16 [l7l (181). The vortic 
ity balance gives trivially VI 2 = 



13 

and is therefore 



omitted for brevity. It remains the following: 



= P-V-B-e 
= V-L-P. 



P, 



(26) 
(27) 



solution via relaxation to the stationary state. 



We neglect the linear molecular friction _L, rewrite (|27|) 
to get P 
following: 

P^B + e 



VP; wie insert this into (26) and get the 

(28) 
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Remarkably, ^> cancels out of (26 27) so that finally 



the oceanographer's standard balance formulation (28) 
is obtained. We come back to this point later. 

In the past experiments with stratified wind tunnels 
gave deep insights into the nature of the turbulence- 
wave interactions (Rohr et al., 1988; Van Atta, 1999; 
Baumert and Peters 2004, 2005). However, they do not 



support relation ( 28 1 . This is the contradiction we men- 



tioned in the Introduction and will be discussed next. 



4.3. The stratified wind tunnel, II 



0, R g > 



4.3.1. The system. This important case is an ex- 
ample of strong horizontal mhomogeneity and non- 
equilibrium conditions. In the entrance facility of the 
tunnel the forced flow passes a fine grid which leaves 
an initial high-frequency short-wave turbulence and 
internal-wave signature in a locally homogeneous fluid 
body of limited size. We consider this fluid body as 
moving with the mean flow in a plug-flow sense. The 
small-scale properties then evolve during the trip within 
the fluid body along the homogeneously sheared and 
stratified tunnel until its end, where the body leaves 
the tunnel, including its small-scale properties. Typi- 
cally an exponential evolution of TKE is observed along 
the longitudinal axis, either exponential growth or ex- 
ponential decay (Van Atta, 1999). The waves were not 
recorded but it must be hypothesized that they did 
reach saturation level. 

This situation is artificial. In a natural hydrody- 
namic system the growth of turbulence is limited be- 
cause it would somewhere begin to reduce its own source 
(shear) by mixing the mean flow (see Fig. [2]) so that an 
equilibrium will sooner or later be reached, correspond- 
ing to the large-scale energy input to the flow. But this 
feedback from the turbulence- wave system to the mean- 
flow system takes time and the travel time through the 
tunnel is too short. 

Exactly this sort of decoupling of processes is what 
the tunnel experimentalists aim at. They wish in par- 
ticular to cut off the feedback P from shear-generated 
wave-energy into TKE, as it contaminates the clear and 
simple picture. In other words, they want to see the 
naked interrelations in stratified shear turbulence as 
studied already by Richardson (1920), Howard (1961) 
and Miles (1961). Those three authors neglected shear- 
generated waves and their subsequent feedback as a re- 
sult of spectral saturation. 

4.3.2. Advection-dispersion-reaction (ADR) 
and the plug-flow concept. At a first glance the 
sheared flow of a stratified wind tunnel seems to repre- 



sent a major problem for detailed analyses. However, 
the Taylor- Aris theory of shear dispersion (Taylor, 1953; 
Aris, 1956; Baumert, 1973; Fischer et al., 1979) allows 
to compute an effective longitudinal dispersion coeffi- 
cient, Dl, and to cast the transport equations corre- 



sponding to (16 17 18) in the following general form 



of an advection-dispersion-reaction equation: 



dY d { dY 



-Y/f. 



(29) 



Here Y is a placeholder for the variables O, /C, and £ , 
and t is an effective time constant of a hypothetically 
decaying (or growing, when f < 0) variable Y. 

As long as the Peclet number of the problem, Pe = 
\(U) 2 t / (A D l)\ ^> 1, it can be shown (e.g. Baumert, 
1973) that the stationary form of ((29j with dY/dt = 0) 
can be simplified into a so-called plug-flow description: 



(U) 



dY 
dx 



-Y/f. 



(30) 



In a wind tunnel this concept is a useful approximation 
because the velocity (U) and thus the above similarity 
number are typically high enough. In a stationary plug- 
flow sense we thus have : 



(U) 

(U) 
(U) 



on 

dx 

WC 

dx 
d£_ 

dx 



1 fS 



, (31) 



= P - $ -B-e + P, (32) 



= * - L - P. (33) 



Now we introduce the dimensionless travel-time coor- 
dinate along the wind-tunnel axis, 

i=x-S/(U). (34) 

We further introduce dimensionless variables via 



Cl = QV2/S, 
JC = K/JCq , 
£ = £/£ . 



(35) 
(36) 
(37) 



Here JCq an d £ o are the initial conditions so that 
)Ct=o — £t—o — 1- Due to the initiation of the flow by 
the grid the vorticity typically begins with high initial 
values, ^t=o ^ 1. 

4.3.3. General case. The above conventions al- 
low to rewrite the transport equations ( [31] |32{ [33| with 
some algebra in dimensionless form: 



dft 

di 



1 

2^r 



1 - fi' 



(38) 
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dJC 
di 
d£ 
di 




(39) 



(40) 



Here we neglected the linear molecular friction term 
L = c\ £ and replaced according to ( 19 1 "J with 
is N 2 /a — B. For B we used {nt and ^W. We now 




1 10 100 1000 



Figure 5. Wind-tunnel turbulence- waves model in di- 
mensionless variables, fl (solid blue) converges soon 
{t « 10) to its structural- equilibrium value, VL^ = 1. K 
(solid red) goes through a minimum and starts at t w 10 
a phase of exponential growth. £ (solid green) rests long 
time close to its initial condition £ o = 1 and enters at 
i k 10 into a phase of exponential growth. The ratio 
P/* (dashed blue, < P < 1) remains initially very 
small but jumps then around t ~ 1500 to its asymp- 
totic value (P/^oo = 1. In this example R g = 0.16, 
a = 10 -6 , and (3 — 1. Notice the double- logarithmic 
character of the presentation. 

abbreviate a — £q/S, /3 = Sq/ICq and a — a x c-i so 
that 



dn 

di 
dt 
di 

dt 
di 



1 


(i-a>; 




2^ 










Ci 2 


M 


7T(T 9 


t 







- a/3£ 2 



(41) 



+ a/3£ 2 , (42) 



(43) 



According to Gregg (1989), in the ocean we have C2 = 
7.4 x 10~ 5 s m~ 2 . a and C2 are the only parameters 



or variables in (41 - 43) which are not dimensionlcss. 
But the two appear only in form of the product a = 
a x C2, which is again dimensionless and was already 
used in ( 42 ) and ( 43 1 . To make an order of magnitude 



3 )- 



2 . This gives 
This value was 



S GM = 0.0036 s"\ and £ w 10" 4 
a characteristic guess of a « O(10" 
used to compute the data in Fig. [5j 

The system (41 - 43) can be solved numerically 



by standard methods if its high stiffness is adequately 
taken into account]^] The first phase of the evolution 
based on ( 41 - 43 ) is shown in Fig. [5]where we may iden- 
tify three regimes with two separating breakpoints. The 
first regime is the initialization or spin-up regime. The 
first breakpoint labels its end at i « 10 and is associated 
with the transition into structural equilibrium. The 
second regime may be called the typical wind-tunnel 
regime (an artificially decoupled or naked turbulence- 
wave system in exponential growth, P w 0). Its end 
is labeled by the second breakpoint at i rj 1500. The 
third regime can be called a hyper-equilibrium regime 
characterized by P = $ and exponential growth, too. 
Both exponential growth regimes exhibit the property 



1 dK 

K di 



constant , 



= 




(44) 
(45) 



We use for regime 3 the name hyper equilibrium because 
it has the condition P = ^> in common with the natural 
equilibrium (oceanographer's regime) but differs with 
respect to the dynamic state: while in the natural case 
all first time derivatives vanish, in the hyper case they 
are constants > 0, but the second derivatives vanish. 

The simulations show that in an ideal stratified wind 
tunnel of infinite length the hyper-equilibrium regime 
can in principle be reached, e.g. for t 3> 1500. But 
such a length can hardly be achieved in practice. Thus 
the most interesting and scientifically unique part in a 
stratified wind tunnel (or a salt-stratified channel flow) 
is the section between the lower and the upper break- 
point. 

In a somewhat sloppy form we may summarize: the 
total travel time through the wind tunnel is too short 
for the shear-generated waves to develop a saturated 
spectrum which would almost equate the breaking term 
P with the generation term 



9 Numcrical overflows may occur. If automatic stiffness tech- 
niques are applied the solution may begin to switch periodically 
between different methods. It is always helpful to reduce the 
maximum time step as far as possible. 



estimate of a we take for 5 the Garret-Munk value, 
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5. Stably stratified natural shear flows: 

R g > o, n > o, s 2 > o 

We now come back to the system (Jl6j [l7j [18] ) and study 
its the TKE balance where the waves and the vorticity 
are assumed to stay in a steady state: 



o = s 2 /2 - n 2 - 

o = n + v-L 

dK 

~dt 



S 2 , 
-P, 



P + P- ^- B- e. 
5.1. TKE balance and mixing efficiencies 



(46) 
(47) 

(48) 



The definition (21 ) allows now to rewrite the TKE bal- 



ance (|48|) as follows, 

(i-/)n + p 



dJC 
~dt 



B 



(49) 



Shear and occasionally overturning waves are qualita- 
tively different generation mechanisms with differing 
spectral signatures, differing mixing efficiencies and dif- 
fering buoyancy fluxes. Both mechanisms are associated 
with buoyancy fluxes which add up (parallel circuitry??) 
to the total buoyancy flux B: 



B 



B' + B. 



(50) 



The so-called mixing efficiency, T, is generally defined 
as T = B/e. Due to the additive nature of B = B' + B 
also r is additive: 



£ 



B' + B 



B' 

e 



B =T> 

£ 



(51) 



With these preliminaries we may rewrite (49 1 as follows: 

dK. 



dt 



(1 - /) II + P - B' - B - e - /* . 



(52) 



5.1.1. Mixing efficiency of purely shear- 
generated turbulence: V = T'(R g ); B = 0, P = 
0, II = 0, S 2 = 0. Using exclusively definitions like 



( 10 1 and rt8| given in the previous sections and subsec- 



tions, we begin here to specify B' as follows: 



B' = £ T' = u N 2 = - N 2 = — 
a a 



1 - 



(53) 



We solve ( 53 ) for T' and get the following 

r' 



a £ 



\ 

00 



f 



(54) 



In the pure shear-generated case, II = S = 0, we thus 
have fl 2 = S 2 /2, and with R g = N 2 /S 2 ([54} rewrites 
as follows: 

r' = — R g {l-2R g ). (55) 

This function is presented in Fig.[6]and is well supported 
by observations as shown in an earlier study by Baumert 
and Peters (2004, 2005). 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 



Figure 6. The mixing efficiency V according to equa- 



tion ( 55 ) for shear-generated turbulence. 



5.1.2. Mixing efficiency of purely wave- 
generated turbulence: r = To = 0.2; P = = 
0, S = 0. The mixing situation is qualitatively differ- 
ent in the case without mean-flow shear where turbu- 
lence is exclusively generated by wave spectra fed even- 
tually by large-scale long- wave external sources. Ac- 
cording to the state of the art (Osborne, 1980; Oakey, 
1982; Gregg et al., 1986; Peters et al. 1988), the mixing 
efficiency may still be taken to be a universal constant, 
r = Tq = 0.2. Clearly, in purely shear-generated tur- 
bulence a gradient-Richardson number is not available 
so that r cannot be a function of it. Further, in this 
case this number is not even defined because the shear 
is zero. 

5.2. The remote-wave closure and 
oceanographer's balance: 

p = * = o, s = o, n > o, s 2 > o 



Here we will derive an approximate relation between 
the large-scale long- wave source of wave-generated tur- 
bulence ("remote waves"), II, and the correspondi ngly 
induced component S 2 in the vorticity balance ( 46 1 . 
We choose a flow without mean-flow shear, S — and 
B' = 0, such that, according to (Eel, fl 2 = The 



TKE balance ( 52 ) reads in this case 



dJC 
dt 



= (l-/)n-B-e, 



(56) 
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and in the steady state: 

(l-/)n = B + e. 



(57) 



The remaining rate / II describes the linear or molecular 
energy loss of the wave field. 
We replace B with T e and get 



(i-/)n = 7£ = 7^r2 2 = 7J/5 2 



(58) 



so that the following is finally the closure of our prob- 
lem: 



S 



2 = (i-/)n 

7 v 



(59) 



6. Natural coexistence equilibria 

We consider now the general natural coexistence of 
waves and turbulence where < / < 1, II > 0, P > 
etc. 

6.1. Vorticity 



We insert the approximate closure relation (|59|) in the 
vc 
1 



vorticity balance ( 46 ) and get with the abbreviation 7 = 

r s 2 | (i-/)n 

2 



n 2 



7 V 



(60) 



For convenience reasons we introduce r\ and choose for 

(61) 



( 60 1 the following presentation 

n 2 



V 



s 2 



1 X 
2 



~ J 

7 



where for our convenience 



* = (W) = (!-/)?■ 
v b 2 - P 



(62) 



6.2. TKE 

Now we discuss the steady-state version of the general 



TKE balance (52) 



{l-f)n + P = B' + B + e + f$. 



(63) 



According to our previous discussions, it may be written 
as follows: 

(l-f)U + uS 2 = -N 2 +je+^-(^-) e. (64) 

<7 <T \ " / 

We rewrite this equation identically as 



{l-f)U+vS 2 = - N 
a 



1 + 



/ (N 



(7q \0 



vVL 2 , (65) 
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and divide both sides of ( 65 1 by v S to get with ( 62 ) 

4' 



1 + X = ^ 
a 



1 + 



f (N 



n ■ 



(To V ^ 

We now expand (N/Q) 4 into (N / S) 4 / (ft / S) 4 



replace rj with (61) and get finally 

2 



7 + 2X 



(66) 



(67) 



Finally, with the helpf of (61), (67) and with the ab- 



breviation ( 66 1 can be brought into the following final 
form: 



l + X 



a \ (T r)' 



so that 



l+X 



Ra 



in- 



^0 V 



(68) 



(69) 



where we remember that 77 = 1/2 + X/j. We solve (69) 
for cr: 



o-(l) 



Ra 



1-7/2- 47/ JS/(7 + 2X) 



(70) 



This function should also satisfy the definition ([8]) of 
the turbulent Prandtl number. This means that 



a(2) ~ 1 - (N/n) 2 ~ n- Rg 



(71) 



and with (61) we get 

a = ( ^ + 2X)fT ° = a m (X, R g ) . (72) 
^ + 2X-2^Rg (2)V ' 9 ' 17 

Here (To and 7 = 1 + Y = 1.2 are a universal constants. 

Now our unknown X = (1 — f)H/P is easily deter- 
mined by equating ( 69 1 and ( 72 ) : 



a {1) (f,X,R g )=a {2) (X,R g ). (73) 
This gives the solution 

X = X(f,R g ). (74) 
Now the knowledge of X as function of / and R g allows 



to present ( 72|) in the following form: 
a 



a (l){ X (f,R g )iRg) = <7 (2)(f,Rg) ■ 



(75) 



The last expression is a family of curves giving us, in 
the way we are used to, for each value of f £ (0, 1) one 
curve a = <j(R g ) as function of the gradient Richardson 
number. 
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6.3. Validity limits 

For simplicity we analyze the system's behavior for / = 
and solve the following equation, 



a (1) (f = 0,X,R g )=a (2) (X,R g ), (76) 



and find 



2R g = ( 7 + 2X)a Q 
2-7 7 + 2(X - ^Rg) 

which is easily solved for X: 



(77) 



X = - 



1 R g -2^R 2 g -a (2-7)7 
2" a (2-j)/2-R g 



Fig. [7] illustrates the function X = X(R g ) around the 



singular point, R g = R 



(i). 



g — -tig 
-^0 



4 1) = (2- 7 )y = (l-f)y=0.2. (79) 



We see that 



X = 



< if R g <R { g\ 
> if R g >R g 1) . 



(80) 



We remember that according to (62 1 negative X <xU/P 



such that X < mean sucking internal-wave energy 
out of the wave-energy pool. We accept therefore that 
for R g < Rg 1 ^ a physically reasonable solution for an 
equilibrium coexistence of waves and turbulence does 
not exist. 

The function X in Fig. [7] exhibits an obvious mini- 
mum. We take (78), differentiate and find the zero of 
dX/dR g here: 



R ( n ] = (2 - 7) o-o = (1 - f ) (7 = 0.4. 



(81) 



Considering this point in Fig. [7J the solution on the left 
is physically unrealistic because a decreasing R g would 
lead to an increase in X. With decreasing R g into a 
region left of the minimum in X we approach a zone 
where the model is simply no longer correct. The reason 
is surely the vorticity balance which has been deduced 
from wind-tunnel experiments and then combined with 
the new remote- wave closure (59). Consequently its 



work is guaranteed only in two cases: 

• purely shear-generated turbulence with exclusion 
of eigen-wave feedback, i.e. P = (wind tunnel); 

• shear-generated turbulence with inclusion of 
eigen-wave feedback, P = (1 - /)(n + > 0, 
but dominant presence of remote waves, i.e. X oc 
TL/P ^> 1 (ocean, atmosphere). 



Unfortunately, the important case with of eigen-wave 
feedback, P = (1 — /)W > 0, but without remote wave- 
energy source, X oc TL/P = 0, is not understood so far. 
It plays a role when a flow setup is shielded from outer 
influences. This may play a role in technical systems 
like circulating cooling ponds where the waves stem not 
from geophysical sources and where they evolve into a 
sufficiently saturated spectrum. This will be discussed 
in the next Section 

6.4. Asymptotic of a for R g — > 00. 



(78) We take the expression ( 78 1 for X and look at very high 



R g where in R g quadratic term dominate constants and 
linear terms such that 



lim X(R g 



-2R a 



1R9 



(82) 



The study of the Prandtl number function at high R q 



is still more easy. We take the left-hand side of (77), 
07]\ = 2 Rg/ (2 — 7), and use the definition 7 = 1 + T = 
1.2 given already previously. Obviously we have 

2 

lim o-(R g ) = =• R g = 2.5 R g . (83) 

h 9 ->oo 1 — r 

For comparison with the numerical solution and with 
observational data we refer to Fig. [9] 




Ra 



Figure 7. The relative wave-energy input to the TKE 
pool, X = (1 — f)U/P, as a function of the gradi- 
ent Richardson number, R g . Negative X indicate with- 
drawal of energy from the wave-energy pool and are 
physically irrelevant. 



7. Geophysically shielded systems 

These systems have no external sources of internal- wave 
energy but they are in coexistence equilibrium of tur- 
bulence and saturated waves with significant eigen-wave 
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feedback, /< 1. Here II — and the steady-state TKE 
balance reads in our notation as follows: 



P = B' + B + e+ fty 

The vorticity balance is 

e_ _ O 2 _ 
P ~ S 2 ~ 



B' + je + /#. 



(1 



1 



(84) 



(85) 



where for young waves, / w 1, the effect of Y is negli- 
gible and the situation is close to the wind-tunnel case. 
Not so for older wave spectra. 



We combine ( 84 1 and ( 85 1 and get after some algebra 



the following steady-state condition, 

2 = ^(l- ¥T |^ 7J y)+7+2(l-/)y, (86) 

where Y appears as a function of the steady-state 
Richardson number R s g , and of the 'wave age' 1 — /. 



The solution Y = Y(R s g ) of ([88]) is presented in Fig. M 
for / = 0. 

The above means that with Y we have a tunable 
parameter which allows us to adjust our model value 
for Rg according to measurements or observations. Un- 
fortunately these are rare for shielded conditions de- 
scribed above so that we are inclined to choose ac- 
cording to tradition the value Rg = 1/4, which cor- 
responds to Y = 0.136. Another choice would be the 
minimum value of Y which is 0.135 and corresponds to 
Rg = 0.265. This somewhat arbotrary situation prob- 
ably explains the large scatter in the measurements of 
a and underlines umso mehr the necessity of dedicated 
experiments and observational programmes. 

We note in closing this Section that the more general 
form of the steady-state vorticity balance is 

e fi 2 _ 
S 2 ~ 



V = n = 



X + (l-f)Y 



1 



(87) 



where the external forces X <x H/P appear together 
with the feedback via eigen waves, Y. Our values of Y 
are situated well below the validity limit of X. In the 
geophysically shielded case (X — 0) and young waves 
(/ sa 1) the action of Y is screened and we have r\ = 
1/2, which is the wind tunnel situation. In the same 
case with adult spectra (/ <C 1) we have the classical 
shielded case. 

8. Application 

In the previous Sections of this report we have looked 
at special physical situations which we knew and under- 
stood sufficiently well. Now we put these pieces together 




0.0 0.1 0.2 0.3 0.4 0.5 



Figure 8. The parameter Y as a function of the 
steady-state gradient Richardson number R g l and for 
/ = 0. The horizontal and the right vertical dashed 
gray lines cross at the minimum of the function Y. 
The left vertical dashed line labels R g l = Y4. Values 
Y ~ 0.135 . . . 0.15 would correspond to mathematically 
0.15... 0.3. 



admissible values Rf 



and write down the full equation system for applications 
also to unknown situations. We choose here the philos- 
ophy of non-stationary embedding. At a first glance 
the associated evolution equations ( 88 - 99 ) look volu- 



minous compared with the lean system of their purely 
algebraic steady-state counterparts, but non-stationary 
embedding avoids stiffness problems right on the most 
fundamental level and is thus substantially more robust 
in the computational practice. 

8.1. Generalized equations 

For an effective notation we define the following differ- 
ential operator, 



V = 

V = 



d 
dt 
d_ 
dt 



d d_ 
dz dz 
d_ . d_ 
dz dz 



(88) 
(89) 



so that the general set of balances can be written for a 
stratified water column as follows: 



vn 


1 


'S 2 


7T 


2 


VIC 


= n 




V£ 


= n 




fit) 






ci£ + 



1-/ 



n 

p + 

a 

C'2 



\ ) s 2 -n 2 
(1 + r)n 2 



C-2 



£2 ' 



,(90) 

(91) 
(92) 
(93) 
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a 


i - n 2 /n 2 ' 


(94) 


P 


= vS 2 , 


(95) 




P S" 2 

R 2 

a n 2 9 ' 


(96) 


V 


= /C/ttQ, 


(97) 


/' 


= v/cr, 


(98) 


R 9 


= N 2 /S 2 . 


(99) 



thus not automtically applicable to atmospheric condi- 
tions. For the latter a comparable result is unknown. 

Gregg (1989) presented a summary of comprehen- 
sive, extensive and direct dissipation and stadardized 
shear observations (S w ) made in ocean waters around 
the globe, where the following conditions applied at 
least approximately: 

• There was almost no mean- flow shear, S ~ 0. 



To apply this theory in form of a numerical model it 
needs to be combined with a scheme which provides us 
with the mean-flow variables from which we may de- 
rive the shear S and the Brunt- Vaisala frequency, N. 
Furthermore we need initial conditions for all variables. 
But the hydrodynamic system is highly dissipative and 
'forgets' the initial conditions soon such that here "rea- 
sonable guesses" would suffice. 

8.2. The parameters 

8.2.1. Overview. T = 0.2 and <tq = 1/2 are uni- 
versal constants. / is a function of time. Its final equi- 
librium value, /oo, depends on the molecular- viscosity 
parameter c\. In many cases it is sufficient to set 

/oo = 0. 

With the exception of i>, the whole system contains 
only one tunable paramter, Y ~ O(10 _1 ). The inner 
physical structure of Y (and T) we have not yet fully 
understood. I.e. we are not able to derive its value from 
other than pragmatic arguments like 'it works', because 
it gives the right gradient Richardson number for the 
geophysically shielded steady-state. 

About the spatial 'diffusivity' of wave packets, v, we 
only know that it scales with the characteristic group 
velocity (c g ) of the wave packets, multiplied by a charac- 
teristic length scale which is possibly the characteristic 
wave lenght of a packet: 



x L. 



(100) 



The last open problem to be discussed is the role of II. 
According to Gregg (1989) it can be estimated from the 
r.m.s. 10-meter wave shear, (S 2 ), as function of time 
and of the region of the world ocean under study. This 
will be done in the next Subsection. 

8.3. Wave-induced dissipation: the ocean case 

The above results apply to stratified oceanic and at- 
mospheric flows as well. The following estimator of 
the long- wave, non-local ("external") energy source II 
is based on extensive studies in the world oceans and is 



• The IGW field was almost perfectly saturated, 
/«0. 

• The observations were done for conditions of 
quasi-steady state, dCl/dt = dfC/dt = d£/dt = 0. 

Gregg established the following empirical relation be- 
tween the wave-induced dissipation rate, e, the 10-meter 
high-pass filtered vertical shear, Sio = A[//10m, and 
the effective Brunt- Vaisala frequency, (TV 2 ) 1 / 2 : 



ai 



(N 2 ) (St ) 



N 2 



c4 

°GM 



(101) 



Very low frequencies have been removed from Sio by 
filtering. The average in (Sf ) is taken over longer ob- 
servation periods. 

Sgm used in ( 101 1 is the so-called Garrett-Munk 
shear: 



a.2 



(N 2 



with the following empirical parameters: 

-10 



Oi = 

a 2 = 
N = 



7 x 10 _iu m 2 s" 3 



1.66 x 10 
5.2 x 10 



- 10 s- 4 . 



~ 3 s" 1 



We insert ( 102 1 in ( 101 1 and get 



«3 



(sio) 

(N 2 ) 



where 



rtl 
« 2 



N£ = 



1.14 x HT 4 n^s" 1 



(102) 



(103) 
(104) 
(105) 



(106) 



(107) 



We further take into account that (see Gregg, 1989) 

(St ) = 2(S 2 W ) 2 (108) 



such that ( 106 ) may be written as follows, 



0,4 



(S 2 10 ) 2 
(N 2 ) 



Ct4 



R a 



(109) 
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with the wave-based gradient Richardson number, 



(S 2 



10/ 



(N 2 



(110) 



The latter needs to be unterschieden von the mean-field 
based gradient Richardson number R g = (N 2 )/(S) 2 . 
The parameter a 4 reads as follows: 



«4 



2a 3 



2.3 x 1CT 4 m 2 ^- 1 



With ( 109 1 we have a solid estimate of II as 
= (l + f)e 



n 



a 4 (l + f) {S2lo) 



(ill) 



(112) 



so that herewith the model system is completely closed 
and we can begin with simulations, provided we have 
a certain knowledge about the wave-field parameters 
(S 2 ) and R g for our ocean or atmosphere region of 
interest. 

9. Comparison with observations 

During the about 15-years walk towards the system 



99 ) partial solutions were step by step con- 



fronted with corresponding observational and experi- 
mental data. The totality of these theory-observation 
comparisons cannot be repeated here. We refer for 
the wind-tunnel experiments to Baumert and Peters 
(2004, 2005), for the Monin-Obukhov boundary layer 
to Baumert (2005a), for the neutrally stratified case to 
Baumert (2005b, 2012). These specific cases are well 
described by subsets of (88 - 99), corresponding indi- 



vidually to the special physical situation studied. 

With respect to a theory-data comparison the 
present article thus concentrates on the case R g ^> 1 
where without mentioning it R g means here always an 
equilibrium or steady-state value. We look namely on 
the relation a = a(R g ), analyzed and summarized by 
Zilitinkevich et al. (2008). These authors used more or 
less the same comprehensive data base for their discus- 
sions like Galperin et al. (2007), Canuto et al. (2008) 
and Kantha and Carniel (2009). In the core these are 
the CASES-99 (stable nocturnal BL, see Poulos et al., 
2002) and the SHEBA experiments (arctic BL over ice; 
see Grachev et al., 2005, 2006). 

Zilitinkevich et al. (2008) conclude that the class- 
wise average of the data is best described by a = 
0.8+5 R g shown in full red in Fig. [9] In our view this line 
is somewhat above the CASES-99 data and our alter- 
native theoretically derived relation, a = (1 + 5R g )/2, 
fits better. It is given in full green in Fig. [9] and located 



somewhat below the full-red curve. Of course, in view 
of the enormous scatter of the original data, the red and 
the green lines are well within the huge overall scatter 
range. 

The two blue lines (full and dashed blue) in Fig. 
[9] represent our simulation results in good agreement 
with the observations of CASES-99. At the same time 
they might illustrate why natural data exhibit such a 
strong scatter. While the dashed blue line is the case 
/ = 0.5 %, the full blue line stands for / = 0. I.e. 
very small absolute variations in / cause large varia- 
tions in the solutions. Generally the huge scatter might 
be caused by the varying age (or degree of saturation) of 
the waves involved. These waves are mostly of non-local 
origin so that in principle all points of the distant space 
are candidates for their birth places and any age of ar- 
riving waves can be expected at our study site. On the 
one hand structural equilibrium needs to be achieved to 
avoid such a scatter. This appears to be difficult un- 
der the action of waves which also modify the flow. But 
still more important seems to be the different age of the 
incoming wave spectra. 

According to our theory Fig. [9] shows also the two 
new critical gradient Richardson numbers for the wave- 
turbulence coexistence in form of two vertical thin black 
dashed lines at R g 1] = (1 - f)/4 = 0.2 and R ( g 2) = 
(1 - f)/2 = 0.4. The higher value labels the lower 
applicability limit of our theory. 

10. Discussion 

This report may be seen as an experiment to draw 
a (dotted) theoretico-physical bottom line under more 
than 50 years of ambitious observational programs and 
experimental laboratory research into stratified small- 
scale geophysical and engineering turbulence and mix- 
ing processes. Looking backwards, today these pro- 
grams appear as a single planned initiative wherein mu- 
tually supplementing pieces fit perfectly together. On 
the water side the efforts took place in the United States 
and the United Kingdom. On the atmospheric side also 
multi-national efforts need to be acknowledged. Our 
results may encourage those who believe that the "tur- 
bulence problem" (Hunt, 2011) is eventually solvable. 
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